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The partition function of an oscillator disturbed by a set of electron particle paths has been 
computed by a path integral method which permits to evaluate at any temperature the relevant 
cumulant terms in the series expansion. The time dependent source current peculiar of the semi- 
classical Su-Schrieffer-Heeger model induces large electron-phonon anharmonicities on the phonon 
subsystem. As a main signature of anharmonicity the phonon heat capacity shows a peak whose 
temperature location strongly varies with the strength of the e-ph coupling. High energy oscillators 
, are less sensitive to anharmonic perturbations. 
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H-H I. INTRODUCTION 

' 

There is at present a growing interest in electron- 
phonon non linearities also triggered by the signatures 
of a large anharmonicity observed in metal diborides [1] . 
Phonon anharmonicities have a long story in solid state 
physics closely related to the inelastic neutron scatter- 
q . ing theory [2,3]. As the anharmonic effects are generally 
ZJ ' small in crystals up to room temperature, second order 
perturbation theory suffices to determine the lifetime due 
I ' to three phonons decay processes while the renormaliza- 
tion of the frequencies, together with the three phonons 
(■■<■) ■ terms, also requires computation of the four phonons ver- 
| tex in the first order diagram plus a (dominant) contribu- 
^sO • tion due to the thermal expansion of the crystal [4] . It is 
' known [5,22,7] that most properties of real materials can 
be well described by replacing the anharmonic phonons 
with the temperature dependent renormalized harmonic 
phonons, that is assuming quasi-harmonic vibrational 
^3 models [8]. Instead, the damping of some anomalous 
bulk and surface phonons requires explicit computation 

1 ■ of the anharmonic interactions [9,10]. However, neither 
'"O ' first principles calculations of anharmonicities based on 

C . density functional perturbation theory [11] nor empirical 
force constant approaches do separate the bare phonon- 
phonon interactions from the non linearities due to the 
electron-phonon coupling. The latter contribution is gen- 
erally incorporated in the former by fitting the third and 
fourth derivatives of the interatomic potential to experi- 
mental thcrmoelastic properties. In systems such as poly- 
mers, whose thermal and conducting behavior is shaped 
by the strength of the e-ph coupling, one would desire to 
estimate the amount of e-ph anharmonicity which may 
become relevant in the low temperature region where, 
instead, the phonon-phonon interactions tend to vanish. 

Theoretical investigations on polymers usually depart 
from the Su-Schrieffer-Heeger (SSH) Hamiltonian [12-14] 
in which the e-ph coupling is the derivative of the elec- 
tron hopping integral with respect to the intersite atomic 
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displacement. As the electron propagator couples to the 
oscillator displacement the SSH model can be attacked 
by the path integral method [15-17] which allows one to 
derive the time dependent probability amplitude for a 
particle in a bath of oscillators [18]. On the other hand, 
considering the electron particle path as the disturbing 
source for the phonon subsystem one may evaluate the 
amount of e-ph anharmonicity by expanding perturba- 
tively the phonon partition function. Analysis of the cu- 
mulant terms versus temperature would permit to assess 
the relevance of the e-ph interactions in the SSH model. 
This paper addresses precisely this issue focussing on the 
computation of some equilibrium thermodynamical prop- 
erties. The path integral method for the model Hamilto- 
nian is outlined in Section II while the results are given 
in Section III. Section IV contains some final remarks. 



II. HAMILTONIAN MODEL 

In one dimension the SSH Hamiltonian is: 

H = ^ Jr,r+l(f}fr+l + fr+lfr) 
r 

J, % r+1 = - 7^ [J + a ( u r - U r+1 )] (0.1) 

where J is the nearest neighbors hopping integral for 
an undistorted chain, a is the electron-phonon coupling, 
u r is the displacement of the atomic group on the r— 
lattice site along the molecular axis, f}, and f r create and 
destroy electrons (i.e., ir band electrons in polyacetylene) 
on the r— group. The non-interacting Hamiltonian is 
given by a set of independent oscillators. We transform 
the real space Hamiltonian of (1) into a time dependent 
Hamiltonian H{t,t') by introducing x{t) and y(r') as 
the electron coordinates at the r and r + 1 lattice sites, 
respectively. The spacial e-ph correlations contained in 
(1) are then mapped onto the time axis by changing: 
u r — > u(t) and u r +i ~~ * U ( T ')- Accordingly we get: 



H 



(r,r') = J T , T '(f\x(r))f(y(r'))+p(y(r'))f(x(T)) 

Jt,t>= -\[J + a{u{T)-u{r'))\ (0.2) 

Eq.(2) shows the semiclassical nature of the model in 
which quantum mechanical degrees of freedom interact 
with the classical variables u(t). As the electron hops 
are not constrained to first neighbors sites along the chain 
H(t,t') is more general than the real space SSH Hamil- 
tonian in (1). Setting r' = 0, u(0) = y(0) = 0, averaging 
the electron operators over the ground state we obtain 
the average energy per lattice site linearly depending on 
the displacements: 

<H N ]> ' = V(x(r))+u(T)j(r) 

j(r) = -a(G[-x{r), -r] + G[x{r), r]) (0.3) 



2 



where N = L/a, with L the chain length and a the 
lattice constant. V(x(t)) (proportional to J) is the ef- 
fective term accounting for the r dependent electronic 
hopping while j(r) is the external source current for the 
oscillator field, G[x(t),t] being the electron propagator. 

Taking a large number of oscillators (ui(r),i = 1..N) 
as the bath for the quantum mechanical particle whose 
coordinate is x(t), the general electron path integral is 
given by: 

<x((3)\x(0) >= Yl J D Ui (r) J Dx{t) 

J\ t (^x 2 {t) + V(x{t)) -£>(t)j(t) 

(0.4) 



xexp 



xexp 



j3 is the inverse temperature, m is the electron mass 
and uoi are the oscillators frequencies. The oscillator 
masses are considered as independent of i, Mi = M and 
hereafter we set M = 10 4 m. After integrating out the os- 
cillators coordinates over the paths Dui(r) and imposing 

a closure condition (x{f3) = x(0)j on the particle paths, 

we obtain the total partition function in the functional 
form: 

Z(j(r)) = Z ph j Dx{T)exp[-^x 2 {r) - V(x(r)) - A{j{r)) 

z ph = Yl 



L 2smh(hwi/3/2) 

h 2 * 1 



A (j( r )) 4M S hut smh(f l Lu l p/2) 
x I* dTj{r) dr" cosh^flr - r"\ -/?/2))j(r") 



The noncquilibrium quantum statistics of the system 
can be derived via (4) through the closed-time path for- 
malism [19] which permits to evaluate dissipative prop- 
erties due to the phonon bath friction [20]. The thermo- 
dynamical properties of the system follow from (5). 



(0.5) 



III. PHONON ANHARMONICITY 



The electron particle path interacts with each of the 
N oscillators through the coupling a (assumed indepen- 
dent of i) of the SSH Hamiltonian. Then, the x(t)— 
perturbed phonon partition function can be expanded in 
anharmonic series as: 
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Z ph [x(r)\ ~ Z ph (l + J2(-1) 1 < C l > x ) 
^ 1=1 ' 

<c k > x = z^ufDu^ji jy Tl 



xexp 



Ui(n)j(Ti) 



(0.6) 



The total phonon partition function Z ph is obtained 
by integrating (6) over a set of perturbing electron paths. 
Since the oscillators are decoupled (and anharmonic ef- 
fects mediated by the electron particle path are here ne- 
glected) we can study the behavior of the cumulant terms 
< C k > x by selecting a single oscillator having energy ui 
and displacement u(t). 

As a first step let's replace the source current of the 
SSH model by a linear (in the electron path) current 
3x{t) = —ax(r). This allows us to derive analytical ex- 
pression for the cumulants as shown, for the lowest order, 
in the Appendix. Let's approximate the electron path by 
its t averaged value: < x(t) >= jj J Q drx(r) = x /a. 
and expand the oscillator path in Np Fourier compo- 
nents: 

Nf 

u{t) = «o + ^] 2^3?m„ cos(w„t) - ^sx n sm(uj n T)j 

71=1 

uj n = 2im//3 (0.7) 
Taking the measure of integration 

„ / 1 \2N F fan ■ ^NfTt) roc 



Np 

X 

n=l 



Nf />oo f-oc 

[I / d$lu n / d%u n (0.8) 

i J —oo J —oo 



with Am = \J irh 2 /3/M, we obtain for the k — th cu- 
mulant the expression 

, r H^ _ rl M^Mjl 1f (2r^) 2 



pfc jfe!7r fc / 2 (w£) fc+1 ± = ± (2mr) 2 + (w/3) 2 
= ax /a (0.9) 

Let's set xo/a= 0.1 in the following calculations thus 
reducing the effective coupling an by one order of mag- 
nitude with respect to the bare value. The trend shown 
by the results hereafter presented docs not depend how- 
ever on this choice since xq/ci and a can be varied in- 
dependently. As the cumulants should be stable against 
the number of Fourier components in the oscillator path 
expansion, using (9) we set the minimum Np through 
the condition 2Npn 3> u>/3. Noticing that odd— k terms 
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vanish the thermodynamics of the anharmonic oscillator 
can be computed by the cumulant corrections to the har- 
monic phonon free energy: 



F<*> = --In 



1 i+E <c ^>J (°- 10 ) 
1=1 



To proceed I need a criterion to find the temperature 
dependent cutoff k* in the cumulant series. In the low 
T limit the third law of thermodynamics may offer the 
suitable constraint to determine k* . Then, given a and lo, 
the program searches for the cumulant order such that 
the heat capacity and the entropy tend to zero in the 
zero temperature limit. At any finite temperature T, the 
constant volume heat capacity is computed as 

c {k) ^ = _ r F ( fe) ^ T + _ 2F W (T + A) + (T) 



f— — 

\A + A 2 



(0.11) 



A being the incremental step and k* is determined 
as the minimum value for which the heat capacity con- 
verges with an accuracy of 10 -4 . Figs. 1(a) and 1(b) 
show phonon heat capacity and free energy respectively 
in the case of a low energy oscillator for an intermediate 
value of e-ph coupling. Harmonic functions, anharmonic 
functions with second order cumulant and anharmonic 
functions with k* corrections are reported on in each fig- 
ure. The second order cumulant is clearly inadequate to 
account for the low temperature trend yielding a nega- 
tive phonon heat capacity below ~ 40K while at high 
T the second order cumulant contribution tends to van- 
ish. Instead, the inclusion of k* terms in eqs.(10) and 
(11) leads to the correct zero temperature limit although 
there is no visible anharmonic effect on the phonon heat 
capacity throughout the whole temperature range being 
Cy ^ perfectly superimposed to the harmonic C' v . Note 
in Fig. 1(b) that the k* corrections simply shift down- 
wards the free energy without changing its slope versus 
temperature. By increasing cxr, the low T range with 

(2) 

wrong (negative) C v broadens (see Fig. 2(a)) whereas 
the k* contributions permit to fulfill the zero tempera- 
ture constraint and substantially lower the phonon free 
energy (see Fig. 2(b)). Thus, for the particular choice of 
constant (in t) source current we find that the e-ph an- 
harmonicity renormalizes the phonon partition function 
although no change occurs in the thermodynamical be- 
havior of the free energy derivatives. Anharmonicity is 
essential to stabilize the system but it leaves no trace in 
the heat capacity. Incidentally we note that a similar role 
for the anharmonic interactions has been proposed to ex- 
plain the low frequency behavior of glasses [21]. Figure 
3(a) displays the k* temperature dependence for three 
choices of e-ph coupling in the case of a low energy oscil- 
lator: while, at high T, the number of required cumulants 
ranges between six and ten according to the coupling, 
k* strongly grows at low temperatures getting the value 
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100 at T = IK for a R = GOmeVA- 1 . The k* versus 
an behavior is pointed out in Fig. 3(b) for three selected 
temperatures. Figure 3(c) shows that by decreasing the 
oscillator energy an increasing number of cumulants has 
to be taken into account to make the thermodynamical 
functions convergent: this is consistent with the expec- 
tations that high energy oscillators are less sensitive to 
e-ph induced anharmonicity. 

Next we turn to the computation of the equilibrium 
thermodynamics of the phonon subsystem in our semi- 
classical SSH model hence, using in (6) the external 
source current given in (3). As the electron propagator 
depends on the bare hopping integral we set J — lOOme^ 
thus assuming a narrow band electron system. Any elec- 
tron path yields in principle a different cumulant con- 
tribution as pointed out in the Appendix. Numerical 
investigation shows however that convergent k— order 
cumulants are achieved by: i) taking Mp = 2 Fourier 
components in the electron path expansion, ii) setting 
for the coefficients {x Ql c TO } a maximum amplitude of or- 
der 0.1 (in units of the lattice constant) and iii) summing 
over <~ 5 2M f+ 1 electron paths. 

As in the previous case, we truncate the cumulant se- 
ries in (6) by invoking the third law of thermodynamics to 
determine the cutoff k* in the low temperature limit and 
by searching numerical convergence on the first and sec- 
ond free energy derivatives at any finite temperature, k* 
does not depend on the specific electron path coefficients 
then, < C k * >= T,( Xo ,c m ) < Ck * >(*„,c m )- Again, we 
can start our analysis from (10) after checking that odd 
k cumulants yield vanishing contributions. Now how- 
ever the picture of the anharmonic effects changes dras- 
tically. The e-ph coupling strongly modifies the shape of 
the heat capacity and free energy plots with respect to 
the harmonic result as it is seen in Figs. 4(a) and 4(b) re- 
spectively. Above a threshold value a ~ lOmeVA -1 (de- 
pending on the energy oscillator here set at u> = 20meV) 
the heat capacity versus temperature curves show a peak 
whose location along the T axis is reported on in Fig. 4(d) 
as a function of a. By enhancing a the height of the peak 
grows and the bulk of the anharmonic effects on the heat 
capacity is shifted towards lower T. At a ~ 60meV 'A -1 
the crossover temperature is around 100K. At high T the 
anharmonic corrections renormalize downwards the free 
energy but their effect on the heat capacity tends to de- 
crease signalling that e-ph nonlinearities are rather to be 
seen in the intermediate to low temperature range. 

Also the k* versus T behavior is much different from 
the previous case (discussed in Fig. 3(a)): it is seen in 
Fig. 4(c) that a few cumulant terms suffice at low tem- 
peratures even at large e-ph couplings while k* grows by 
increasing T for any a. Then, k*(T) does not necessarily 
provide a measure of the degree of e-ph anharmonicity 
as revealed by the heat capacity behavior. At room tem- 
perature we find k* = 14 as the maximum value for the 
range of a's here considered with a low energy oscillator. 
However, by taking u = WOmeV, k* drops (as shown in 
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Fig. 4(c)) attaining the value 6 at around room tempera- 
ture for a = QOmeVA" 1 . This would suggest that also 
in the SSH model high energy phonons are expected to 
be less anharmonic. This idea is confirmed by looking at 
Fig. 4(d) where the crossover temperature in the oscillator 
heat capacity curves is reported on, at fixed a, versus u>. 
The heat capacity peak shows up at u < bOmeV and its 
temperature location is shifted downwards by decreas- 
ing the phonon energy whereas, at larger energies, the 
heat capacity has a positive derivative in the whole range 
of temperatures. Finally we note that recent numerical 
studies of a classical one dimensional anharmonic model 
undergoing a Peierls instability [22] also find a specific 
heat peak as a signature of anharmonicity. 

IV. CONCLUSIONS 

We have studied the equilibrium thermodynamics of 
an electron-phonon system looking in particular at the 
anharmonic effects produced by the electronic subsys- 
tem on the phonon oscillators. The path integral for- 
malism permits to analyse the e-ph non linearities as a 
function of the source current (peculiar of the Hamilto- 
nian model) which disturbs the harmonic oscillator. The 
phonon partition function has been expanded in fc— or- 
der cumulant series and, at any temperature, we have 
determined the cutoff k* that makes first and second free 
energy derivatives convergent. In the zero temperature 
limit the constraint imposed by the third law of thermo- 
dynamics has been fulfilled. Rather than being a unique 
measure of e-ph anharmonicity k*{T) turns out to be a 
model dependent function whose values may vary consid- 
erably according to the physical quantities one decides to 
sample. As a general trend we find however that, at fixed 
T and e-ph coupling, higher energy oscillators are more 
stable and require a lower k* for the computation of their 
thermodynamical properties. 

We have first considered the case of a source current 
proportional to the time averaged electron path and ana- 
lytically derived the cumulant expansion for the phonon 
partition function. This current induces a slope preserv- 
ing shift in the phonon free energy versus temperature 
hence, the heat capacity does not show any correction 
with respect to the harmonic result in spite of the high 
number of anharmonic terms which appear in the cumu- 
lant series. 

Turning to the source current peculiar of the time 
dependent, one dimensional, semiclassical Su-Schricffcr- 
Hecger model we find a striking evidence of e-ph anhar- 
monicity in the constant volume heat capacity of a single 
oscillator. As a main feature the heat capacity exhibits a 
peak whose height and location on the T axis varies with 
the strength of the e-ph coupling and the energy of the os- 
cillator: while strong couplings shift the main body of the 
anharmonic effects towards low T, high energy phonons 
prove to be less affected by e-ph corrections. 
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FIG. 1. (a) Phonon heat capacity and (b) Phonon free 
energy calculated in i) the harmonic model, ii) anharmonic 
model with second order cumulant, iii) anharmonic model 
with k* cumulants (see text), as is the effective e-ph cou- 
pling in units meVA~ and ui is the phonon energy. 



FIG. 2. As in Fig.l but with larger e-ph coupling. 



FIG. 3. (a) Number of cumulants required to obtain a con- 
vergent phonon heat capacity at any temperature for different 
choices of e-ph couplings, (b) Number of cumulants yielding 
a convergent phonon heat capacity at any e-ph coupling for 
three selected temperatures, (c) Number of cumulants yield- 
ing a convergent phonon heat capacity at any oscillator energy 
for three selected temperatures. 
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FIG. 4. Anharmonic (a) phonon heat capacity and (b) free 
energy versus temperature for eight values of e-ph coupling. 
The harmonic plots are also reported on. A low energy oscil- 
lator is assumed, (c) Temperature dependence of the cutoff 
k* for: i) the weakest and strongest e-ph couplings given in 
(a) and (b) with u> = 20meV, ii) the strongest coupling case 
with u> — WOmeV. (d) Temperatures at which the heat ca- 
pacity gets the maximum value (see (a)) versus e-ph coupling 
at fixed ui (bottom axes) and versus ui at fixed coupling (top 
axes). 



APPENDIX 



Let's assume a source current j x (r) = —ax(r) as the 
disturbing term for the oscillator field, a being the e- 
ph coupling and x(t) being the electron path. Then the 
k-th order cumulant term contributing to the x(t)— per- 
turbed phonon partition function is given by 

<C k > 

1=1 L 
M, 
~2 



z ^Jlf Du{T) tl I drMn)x(n) 



(0.12) 



We use the oscillator path and the functional measure 
of integration given in the eqs. (7) and (8) of the text re- 
spectively. Then, the electron path expansion in Fourier 
components 

M F 

x(t) = x a +'^2 ( c m cos(w m r) + d m sm(uj m T)j 

m— 1 
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LU m — 27rm/P 
dm — 2-^5 t r m 



(0.13) 



will be truncated at Mp = Np. Without any loss 
of generality we set d m = and derive the following 
expressions for the lowest even order cumulants: 



< C > ( Xo ,c m ) 
Np 



J P h 2 ^ 



+ 



1 



E 



< c 4 > 



2ir ^ (2irm) 2 + (/3w) 2 

m=l v / v / 

!« 4 (/3A M ) 4 



P h 4! /3w 
2ci 



F(Np) 



(x a ,C m ) 

2tt 2 (/3 W ) 2 ^ (27rm) 2 + (^ 
2ci 



3x 4 



E 



+ 



3 
8^2 



E 

m— l,/>m 

Ay 

47T 2 ^ ((27r77j) 2 + (f3uj) 2 y 



2ir 2 (f3u) 4 



2c? 



(27rm) 2 + (/3w) 2 (2tt/) 2 + (/3w) 2 

„4 



+ 



n 



(27rm)^ 



^ (27rm) 2 + (/3^) 2 



(0.14) 



Odd fc— cumulants vanish at any order. Setting a max- 
imum amplitude for the coefficients x and c m and inte- 
grating eqs.(14) over a class of electron paths one finds 
the total cumulant contributions to the partition func- 
tion. 
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